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Motivation: Detailed experimental knowledge of the level structure of light weakly bound nuclei is necessary to 
guide the development of new theoretical approaches that combine nuclear structure with reaction dynamics. 

Purpose: The resonant structure of *B is studied in this work. 

Method: Excitation functions for elastic and inelastic ^Be+p scattering were measured using a ^Be rare isotope 
beam. Excitation energies ranging between 1.6 and 3.4 MeV were investigated. An R-matrix analysis of the 
excitation functions was performed. 

Results: New low-lying resonances at 1.9, 2.5, and 3.3 MeV in *B are reported with spin-parity assignment 0^, 
2"*", and 1"*", respectively. Comparison to the Time Dependent Continuum Shell (TDCSM) model and ab initio 
no-core shell model/resonating-group method (NCSM/RGM) calculations is performed. This work is a more 
detailed analysis of the data first published as a Rapid Communication. [J.P. Mitchell, et al, Phys. Rev. C 82, 
011601(R) (2010)] 

Conclusions: Identification of the 0"*", 2"*", 1^ states that were predicted by some models at relatively low 
energy but never observed experimentally is an important step toward understanding the structure of *B. Their 
identification was aided by having both elastic and inelastic scattering data. Direct comparison of the cross sections 
and phase shifts predicted by the TDCSM and ab initio No Core Shell Model coupled with the resonating group 
method is of particular interest and provides a good test for these theoretical approaches. 



I. INTRODUCTION 

One of the main goals of modern nuclear theory is to 
combine the nuclear reaction models with nuclear struc- 
ture calculations to provide the unified framework that 
allows the calculation of level spectroscopy and reaction 
cross sections starting from the same established nuclear 
Hamiltonian. Several theoretical approaches have been 
suggested to advance this goal. Broadly, two major direc- 
tions can be identified, phenomenological and ab-initio. 
The first one uses the well established shell-model Hamil- 
tonian and couples it with the corresponding reaction 
channels. The recoil corrected continuum shell model 
(RCCSM) [ij andjhe time-dependent continuum shell 
model (TDCSM) [2] are examples of these approaches. 
The second major direction is the attempt to calculate 
the cross section starting from both bare nucleon-nucleon 
forces and three-nucleon forces. One example of this 
approach is the no-core shell model combined with the 
resonating-group method (NCSM/RGM) [3]. The very 
attractive feature of these developments is that the exci- 
tation functions of the resonance reactions, such as elas- 
tic and inelastic nucleon scattering, (p,n) and (p,a) re- 
actions, etc., can, in principle, be calculated and directly 
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compared to the experimental data. This is in addition 
to all known structure data. However, this comparison 
is not as straightforward as it may appear. Because of 
model space truncation, limitations from nucleonic de- 
grees of freedom and numerical complexity, it is natural 
to expect that the nuclear spectrum at the low excita- 
tion energy is reproduced better than the spectrum of 
the higher lying excited states by any model. Therefore, 
it is desirable to verify the theoretical predictions in the 
region of low excitation energy first and weakly bound 
nuclei provide a good test for these models. Here, the 
continuum appears at low energy, thus permitting exam- 
ination of the structure-reaction transition. Moreover, 
because of truncation of the model space, parameters of 
models are adjusted to the well known spectrum of stable 
nuclei, resulting in unsurprisingly reasonable agreement 
with the experimental data for these nuclei. The bet- 
ter test is provided by exotic, weakly bound nuclei. The 
neutron deficient Boron isotope, ^B, is a prime example. 
Its proton separation energy is only 137 keV and all of 
its excited states are in the continuum, as can be seen 
in its level structure in Fig. [T] In addition, this nucleus 
has been a subject of numerous theoretical studies. In 
the recent ab initio NCSM/RGM analysis [3[ of ^B the 
proton-|-''Be elastic scattering phase shifts as well as the 
cross section for the ^Be(p,p') and the ^Be(p,7) reac- 
tions were calculated. Direct comparison of the exper- 
imental results on the ^Be(p,p) and ^Be(p,p') reactions 
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FIG. 1. (Color online) The level schemes for *B and its mirror 
*Li. States from our previous work [3| are in red. The dashed- 
dotted line indicates that the state is tentative. 



with these calculations and also the analysis of the ex- 
perimental data using the TDCSM approach is the main 
objective of this work. 

The excitation function for ''Be+p has been previ- 
ously measured in several experiments [B-Q- However, 
the thick target inverse kinematics experimental method 
used in all of these measurements did not allow for sepa- 
ration between elastic and inelastic scattering except for 
the data from Q, where measurements were performed 
at energies below the inelastic scattering threshold. In 
Ref. Q, an attempt has been made to use 7-proton co- 
incidence to identify the inelastic scattering events, how- 
ever, the elastic excitation function still appears to be 
contaminated with inelastic events (see section [II] for ad- 
ditional comments) . The experiment described here does 
not suffer from such deficiency because the Intermediate 
Target Thickness Approach has been applied. This ap- 
proach allowed for measurement of a significant fraction 
of the ''Be-l-p excitation function, while simultaneously 
detecting the ''Be recoil in coincidence with protons in 
order to distinguish between elastic and inelastic scat- 
tering events kinematically. Therefore, we did not use 
experimental data from the previous higher energy mea- 
surements d, H, Q in the analysis but we included the 
low energy ''Be(p,p) elastic scattering data between 0.3 
and 0.75 MeV measured in Ref. Q. The subset of the 
data reported here was first published in [1]. 

This paper contains a more detailed description of the 
experimental results and also extends the previously mea- 
sured excitation energy region to higher energies. A de- 
scription of the experimental method that was used to 
measure the excitation functions for ^Be-|-p elastic and 
inelastic scattering between 1.6 to 3.4 MeV is given in 
section [TTl The analysis of the experimental data was 
performed using the multi-channel multi-level R-matrix 
approach and is described in section III. Section IV con- 
tains a discussion of this finding and its consistency with 
the previous experimental data on the ®B and ^Li nu- 




FIG. 2. (Color online) The experimental setup. The ^Be 
beam was delivered by the RESOLUT facility (on the left). 
The C2H4 targets of various thicknesses were used. The pro- 
tons were detected in an array of three Micron Semiconduc- 
tors S2 detectors and the ^Be's were measured in an S2 down- 
stream. (The inset provides a more detailed view of the de- 
tector arrangement.) 



clei and discusses whether it is possible to explain the 
new experimental data without introducing the new res- 
onances in *B. Analysis of the new experimental data in 
the framework of the Time Dependent Continuum Shell 
Model is presented in section V. Detailed comparison 
of the phase shifts extracted from the analysis of the 
p-|-^Be experimental data to the predictions of the ah 
initio NCSM/RGM model is given in section VI. Con- 
clusions are given in Chapter VII. 



II. EXPERIMENT 

The excitation function for p-l-^Be elastic and inelastic 
scattering between 1.6 and 3.4 MeV in c.m.s. was mea- 
sured at the John D. Fox Superconducting Accelerator 
Laboratory at Florida State University. A radioactive 
beam of ^Be was produced using the iH('^Li,'^Be)n re- 
action. A primary ''Li beam was accelerated by a 9MV 
SuperFN Tandem Van de Graaff accelerator followed by a 
LINAC booster. The primary target was a 4 cm long hy- 
drogen gas cell with 2.5 /im Havar entrance and exit win- 
dows. The gas cell was cooled by liquid nitrogen and had 
a gas pressure of 390 mBar. The in-flight production rare 
isotope beam facility RESOLUT was used to separate 
^Be from other reaction products and the primary beam. 
RESOLUT is a set of two superconducting solenoids, 
dipole and quadrupole magnets and a superconducting 
resonator. Three ^Be beam energies were used in this 
experiment: 27.2, 22.0, and 18.5 MeV. The typical in- 
tensity of the ''Be beam was 10^ pps. The composition 
of the beam was 70% ^Be and 30% ^Li contaminant. Di- 
agnostics of the secondary beam were performed using a 
position sensitive Micro-Channel Plate detector installed 
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FIG. 3. Scatter plot of kinematic coincidence between pro- 
tons and ^Be ions. Regions which correspond to elastic and 
inelastic scattering are labeled. 



between the dipole magnet and the second solenoid (2.7 
m before the C2H4 target) and the AE-E telescope con- 
sisting of an ionization chamber (used as AE detector) 
backed by a 50x50 mm^ 16x16 Silicon strip detector po- 
sitioned 66 cm downstream from the secondary target. 

A sketch of the experimental setup is shown in Fig. 
[5J A solid polyethylene (C2H4) target of thickness opti- 
mized for the given beam energy (see description below) 
was used. A set of three annular Micron Semiconductor 
silicon strip detectors (S2 design) for the proton recoils 
were positioned 5, 6 and 7 cm downstream from the tar- 
get, respectively. Another S2 detector for the ''Be recoils 
was positioned 24.5 cm from the target. The S2 detector 
has annular geometry and consists of 16 segments and 
a side of rings that allow for the scattering angle of the 
products to be determined. The first in the set of three 
proton detectors was a AE detector of 65 /im, while the 
other two and the ^Be detector were 500 //m each. 

The target thickness was optimized for maximum en- 
ergy losses of the ^Be ions in the target while ensuring 
that all ''Be recoils make it out of the target with enough 
kinetic energy left to be detected in the downstream S2 
detector. Kinematic coincidence between protons in the 
array of three S2 detectors and the ''Be recoils in the 
downstream S2 detector were then used to identify the 
scattering events. The 65 ^m AE S2 detector was used 
only in the initial stage of the experiment to verify that 
kinematic coincidence between light and heavy recoils are 
enough for clean identification of the p-l-'^Be elastic and 
inelastic scattering events. This detector was then re- 
moved. Measurements at the beam energies of 22 and 
18.5 MeV were performed without the 65 ^m detector, 
while that for the 27.2 MeV energy included the 65 /im 
AE detector. Time between the events in the proton and 
^Be detectors was measured with resolution of about 3 
ns in order to eliminate random coincidence background. 
Elastic and inelastic scattering processes can be distin- 
guished, because complete kinematics of the events are 
measured. More specifically, events that have two equal 



FIG. 4. The excitation function for '^Be-l-p elastic and in- 
elastic scattering at 148±4° and 146±4° degrees respectively. 
Results from runs at three different energies of ■^Be beam are 
shown. The squares correspond to the run at 18.5 MeV of 
^Be with a 1.5 mg/cm^ target, the vertex up triangles are 
data taken at 18.5 MeV with the 2 mg/cm^ target, the cir- 
cles are the 22 MeV data with the 2.5 mg/cm^ target, and 
the vertex down triangles are from the 27.2 MeV run with 
a 2.6 mg/cm^ with solid markers representing ^Be-|-p elastic 
scattering and open markers the inelastic p(^Be,p')^Be(| ) 
scattering excitation functions. 

energy protons would have different energy of ^Be recoils 
if they originate from different (elastic/inelastic) pro- 
cesses. This is due to different reaction Q- value and kine- 
matics, and also effective target thicknesses (and hence 
energy losses) experienced by the heavy recoils. The in- 
elastic events that produce protons with the same kinetic 
energy as elastic events take place earlier (upstream) in 
the target, where a negative reaction Q-value is compen- 
sated by the higher energy of the ^Be projectile. (See also 
Ref. [9|] for details on this experimental technique). The 
2D scatter plot for the kinematic coincidence between 
protons and ^Be is shown in Fig. [31 The kinematic loci 
which correspond to elastic and inelastic scattering pro- 
cesses are labeled and outlined with contours. Polyethy- 
lene target thicknesses used in this experiment were 2.6, 
2.5 and 1.5 mg/cm^ for the 27.2, 22 and 18.5 MeV beam 
energies respectively. In addition, a separate run at 18.5 
MeV of ^Be beam energy was performed with a slightly 
thicker (2 mg/cm^) target, to extend the measured ex- 
citation function to lower energies without changing the 
energy of the beam. Under this condition coincidence be- 
tween the highest energy protons and the ^Be recoils are 
lost (the heavy recoils produced at the beginning of the 
target do not make it through). Only the lower energy 
part of this spectrum was used in the analysis. 

Fig. m shows excitation functions for resonance elas- 
tic and inelastic scattering of ^Be+p measured in four 
different runs. The vertex up triangles correspond to 
the ^Be run at 18.5 MeV with the 2 mg/cm^ target, the 
squares are 18.5 MeV ''Be with the 1.5 mg/cm^ target 
data, the circles are 22 MeV ''Be with 2.5 mg/cm^ data, 
the vertex down triangles are from the 27.2 MeV run 
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with the 2.6 mg/cm^ target, and in all cases, the filled 
markers are for elastic scattering and the hollow mark- 
ers for inelastic scattering. The angular resolution of the 
experimental setup, as determined by the pitch of the 
rings in the S2 detector, distance from the target and 
the size of the beam spot on the secondary target, was 
1.25°. We used binning of 4° in the lab frame, com- 
bining events recorded by 12 rings of the S2 detector 
into one spectrum. Excitation functions at three scat- 
tering angles were obtained this way. These angles are 
148±4°, 140±4° 132±4° in c.m.s. for elastic scattering 
and 146±4°, 138±4° 130±4° for inelastic scattering. Ab- 
solute normalization of the cross section was performed 
using the known excitation functions for ^Li+p clastic 
scattering. These excitation functions were extracted 
from the experimental data using the same procedure as 
for the ''Be-l-p elastic scattering, therefore, by normaliz- 
ing the ''Li-l-p data to the known ^Li-|-p cross section and 
taking into account the ratio of the ^Be ions to the ^Li 
ions in the secondary beam (as measured by the degree 
ionization chamber and silicon strip detector), accurate 
normalization is achieved. Note that this normalization 
procedure automatically takes into account the efficiency 
of the experimental setup. A sample of the ''Li-|-p exci- 
tation function measured in this experiment is shown in 
Fig. [5] (solid circles) and compared to the experimental 
data from [l^, [ll| . Excitation functions extracted from 
our data agree well with the differential cross section for 
elastic and inelastic scattering of ^Be-|-p measured at sev- 
eral energies of ^Be using a thin target approach and 
reported by U. Greife, et aL, [l3|. The excitation func- 
tions of Yamaguchi et al. |8| however, differ from ours, 
especially in the inelastic channel where they found the 
excitation function to be fairly flat across their entire en- 
ergy range measured, while our results have a large peak 
at an excitation energy of 2.5 MeV. This discrepancy may 
be related to the background in the Nal scintillator de- 
tectors used in Ref. jSj] that could have prevented a clean 
7-proton coincidence spectrum to be extracted. 



III. 



R-MATRIX ANALYSIS 



The excitation functions for elastic ^H(''Be,p)^Be(g.s.) 
and inelastic iH('^Be,p')'^Be(l/2-; 0.43 MeV) scattering 
were analyzed using a two channel, multi-level R-matrix 
approach. The natural starting point for the analysis is 
to introduce only the excited states of ®B that were iden- 
tified in previous experiments [Ts], the l"*" at 0.77 MeV, 
the 3+ at 2.32 MeV and the broad 2" at -3 MeV. These 
three states reproduce the excitation function for p+^Be 
elastic scattering between 0.5 and 3.5 MeV reasonably 
well, as shown in Fig. [5] (top) by the solid line. However, 
it is not possible to explain 30 mb/sr inelastic cross sec- 
tion at 2.5 MeV if only known states in ®B are considered 
(Fig. ini (bottom)). This failure can be understood from 
the following simple considerations. The first excited 1^ 
state at 0.77 MeV is too narrow to have any significant 
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FIG. 5. The excitation function for ^Li-fp elastic scattering at 
148±4° is shown with solid circles. This excitation function 
was measured simultaneously with ^Be-|-p (the rare isotope 
beam composition was 70% ^Be and 30% ^Li) and used for 
absolute normalization. The same excitation function from 
[lOl. [Till is shown for comparison with open circles. 



impact on the excitation functions at energies above 1.6 
MeV. The second excited state, 3^ at 2.32 MeV, can 
only decay to the 3/2~ ground state of ^Be because de- 
cay to the l/2~ first excited state requires angular mo- 
mentum of ^ = 3. Therefore, even if the corresponding 
reduced width is large the inelastic partial proton width, 
Fp/ = 2Pf(fci?)7^, would still be small compared to the 
elastic partial proton width due to a small penetrabil- 
ity factor for high angular momentum decay. Hence, the 
cross section for population of the first excited state in 
^Be due to the 3^^ resonance in ^B, determined by the 
TpTpi /FjQ^ ratio, is small. The same is true for the broad 
2" state in *B at «3 MeV as it can only decay to the first 
excited state in ^Be with angular momentum £ = 2 while 
decay to the g.s. proceeds with £ — 0. Fig. |6] shows the 
results of an R-matrix calculation with only previously 
known 1+, 3+ and 2" states at 0.77, 2.32 and 3.7 MeV 
with reduced width parameters evaluated using the TD- 
CSM (more details on TDCSM calculations are given in 
SectionlVj and known total widths of these states. (Exci- 
tation energy and width of the 2~ were adjusted slightly 
to produce a better fit.). It is clear that while the elastic 
scattering data is well reproduced, the inelastic scatter- 
ing data cannot be explained by the known states. 

Based on the level scheme of ^Li (Fig. [T]) it is natu- 
ral to introduce the second 1"*" state in ®B at an excita- 
tion energy around 3 MeV. Reduced widths for this state 
were chosen according to TDCSM calculations carried 
out with the Cohen-Kurath CKI interaction [3]. It was 
verified that these reduced widths reproduce the known 
width of this state in ^Li (-1 MeV). The short dashed 
curve (red) in Fig. [5] shows the effect of the l'^ state on 
the fit. While the elastic excitation function is fitted well, 
the inelastic cross section is still underestimated. Even if 
this state is shifted to 2.5 MeV, where inelastic scatter- 
ing has its maximum cross section, it still underestimates 
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FIG. 6. (Color Online) R-matrix fit of the elastic and inelas- 
tic ^Be+p scattering with known 3'^ and 2~ states, a second 
excited 1^ seen in *Li, and the "background" 1~ state intro- 
duced at higher energy. The solid curve corresponds to only 
3^ and 2~ states at 2.3 and 3.5 MeV, respectively. The red 
short-dashed curve includes the contribution of the higher ly- 
ing l"*" states assumed at 3.0 MeV. Dash-dotted purple curve 
shows the 1+ state shifted to 2.5 MeV and the long dashed 
green curve also includes the 1~ state introduced at 5 MeV. 



the data (dash-dotted (blue) curve in Fig. Finally, in 
an attempt to increase the inelastic cross section with- 
out using new states below 3 MeV we introduced a 1" 
"background" state. This state can decay to the first 
excited state of ^Be with i = 0, therefore it may con- 
tribute significantly to the inelastic cross section. The 
reduced widths for the 1~ state were evaluated using the 
shell model, and the state was introduced at 5 MeV. As 
expected, the 1~ background state increased the inelas- 
tic cross section overall (long-dashed green curve in Fig. 
in])- But even with this state included the inelastic cross 
section cannot be reproduced. 

The ab initio calculations for ^B 0, predict 
three more positive parity (p-shell) states at low exci- 
tation energy. These are the O+i, I+2 and 2+2- The 
excitation energies for these states vary between 2 and 
6 MeV depending on the three-body force parameteriza- 
tion and the specifics of the calculations. Similar results 
are obtained in shell model calculations (excitation en- 
ergies of these "missing" states vary between 2 and 6 
MeV in the shell model as well, depending on the resid- 
ual interaction used). Therefore, it is natural to intro- 
duce these states in an attempt to reproduce the large 
inelastic scattering cross section. The I+2 state has al- 
ready been introduced. That leaves only the O+i and 
2+2- Introduction of a new 2+ state placed at 2.5 MeV, 
reproduces both the magnitude and angular dependence 
of the observed peak in the inelastic cross section while 
keeping the elastic excitation function in agreement with 
the experimental data (green dash-dotted curve in Fig. 
[7]). However, even with this new state the cross section 
for inelastic scattering below 2.3 MeV is still underesti- 



mated. The 2+ state should have a relatively small width 
(270±40 keV) to fit the observed peak-like structure in 
the inelastic excitation function at 2.5 MeV and its in- 
fluence below 2.3 MeV is small. Introducing the 0+ state 
at an excitation energy of 1.9±0.1 MeV with a width of 
5301^00 ^sV allows the inelastic scattering data to be fit 
below 2.3 MeV without destroying the fit to the elastic 
scattering data (solid blue line in Fig. [7]). It was verified 
that a 1+ spin-parity assignment for this state is not pos- 
sible as it will ruin agreement with the elastic scattering 
data (red dashed line in Fig. [3. 

The low-energy data from Q were used (Fig. [SJ to 
provide additional constrain on the behavior of the phase 
shifts at low energy. It proved to be particularly impor- 
tant for the negative parity phase shifts. We used the 
predictions of the ab initio calculations ITJ for the 2~ 
and 1~ phase shifts as the starting point, but the best fit 
was achieved with the negative parity phase shifts differ- 
ent from [Trj . (It is discussed in more detail in Chapter 

f) The best fit that included the low energy data from 
and data from this experiment was achieved using R- 
matrix parameters given in Table HI The normalized 
for the best fit was 0.89. States shown in parenthesis in 
Table U are the broad "background" states that are used 
in R-matrix formalism to produce the desired behavior 
of the corresponding phase shifts. 



IV. NEW STATES IN LIGHT OF PREVIOUS 
EXPERIMENTAL DATA 

Here we study/assess if the low lying 0^ and 2 J states 
are consistent with the available experimental data on 
^B and ^Li nuclei. The structure of ®B has been ex- 
tensively studied in p-|-^Be resonance elastic scattering 
experiments 0411 ■ In Ref. Q the ^Be-|-p excitation func- 
tion of elastic scattering was measured in the cm. energy 
range from 0.3 to 0.75 MeV. The new states are at 1.9, 
2.5, and 3.3 MeV excitation energies and their influence 
on the low energy part of the excitation function is very 
small. In general, the fit to the elastic scattering data 
does not require the low lying 0^ and 2^ states. The ex- 
perimental data in Ref. [5, 6, 8] were fitted with only 1+ 
and 3+ states at 0.77 and 2.32 MeV and a 2~ state at ^ 3 
MeV. (Presence of a 1'^ state at 3 MeV was suggested 
in Ref. 5].) However, the new states have little influence 
on the excitation function for elastic scattering. There- 
fore, the absence of these states in the R-matrix analysis 
of the elastic scattering data cannot be used as an argu- 
ment against these states. It is interesting to note that in 
all three previous measurements 0, H, 0] the cross section 
at the resonance energy of the 3+ state was measured to 
be « 190 mb/sr at 180°. The R-matrix fit to our elastic 
scattering data produces a lower cross section at 180°, 
«160 mb/sr. This is not surprising because the experi- 
mental technique used in previous measurements did not 
separate elastic from inelastic scattering. Protons from 
inelastic scattering were contributing to the "elastic" ex- 
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FIG. 7. (Color online) Elastic and inelastic excitation functions for '^Be+p scattering. The green dashed-dotted curve is a fit 
with the previously known l""", 3'*', and 2~ states, as well as a 2^ at 2.54 MeV to reproduce the peak in the inelastic data. The 
red dashed curve is a fit with the aforementioned states and the state. The blue solid curve is the best fit with all states 
from Table m The black dashed-dotted curve is a TDCSM calculation. 



TABLE I. Parameters of resonances in *B from the R-matrix best fit. States in parenthesis are outside of the measured 
excitation energy range but provide essential "background" through low energy tails. The energy eigen-value and the reduced 
widths amplitudes for ^Be(p,p) and ^Be(p,p')^Be(l/2~) scattering with channel spins 1 and 2 for the former and and 1 for 
the latter used in the R-matrix fit are shown in columns 6-10. We used 4.20 fm as the channel radius for both the elastic S=l,2 
and inelastic S=0,1 channels. 
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FIG. 8. (Color online) Excitation function for ^Be+p elastic 
scattering at low energies from [J at 123.6° and 163.8°. The 
best fit is a solid red curve. The calculated cross section 
was convoluted to account for 30 keV experimental resolution 
reported in Ref. [J. Systematic errors were included into 
the error bars. The dashed green curve is the R-matrix fit 
with the ab initio 1"*" phase shifts from [TJ for the 1^ state 
(channel spins 1 and 2 contribute about equally). 



citation functions which resulted in higher cross section 
values of the measured "elastic" excitation functions. 

More experimental information is available regarding 
the structure of the mirror nucleus, ^Li. One and two 
neutron transfer reactions, ''Li(d,p) [l^ and ^Li(t,p) [l^ 
were used to populate states in *Li. It is very unlikely 
that bound states in *Li could have been missed in these 
experiments. Therefore, the 0^ and 2^ states are prob- 
ably above the neutron decay threshold (2.03 MeV) in 
8Li. 

The excitation function for the ^Li(n,7)*Li reaction 
was measured at low cm. energies (up to 1 MeV) [13- 
Only the 3+ state at 2.25 MeV (0.22 MeV above the 
neutron decay threshold) has been observed. In princi- 
ple, lack of evidence for the 0+ and the 2+ states in the 
'^Li(n,7)^Li excitation function cannot be considered as 
a decisive argument against their presence. If the partial 
7 width (F-y) for these states is small then they can be 
hard to identify within the background from direct neu- 
tron capture and resonance capture due to the 3+ state. 
Fig. [5]shows TDCSM calculations of the (n,7) excitation 
function with the known 3"*" and 1^ states and the new 
0+ and 2+ states at 2.4 and 2.5 McV (top panel) and at 
2.8 and 3.3 MeV (bottom panel). It is clear from this 
figure that observation of the new states in the ''Li(n,7) 
reaction is difficult. 

Resonances in ^Li at excitation energies of up to 9.0 
MeV have been studied in elastic and inelastic n-f^Li 
scattering and analyzed using the R-matrix approach 



Ref. 



24| . where the new low- lying states were sug- 
gested. For example, the 0+ state at 3.02 MeV was intro- 
duced. Unfortunately, the n-|-^Li excitation function is 
relatively featureless, which makes R-matrix analysis am- 
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FIG. 9. (Color online) The ^Li(n,7) reaction excitation func- 
tion calculated using TDCSM approach. The known 3^ and 
1^ states and the new 0^ and 2"*" states at 2.4 and 2.5 MeV 
(top panel) and at 2.8 and 3.3 MeV (bottom panel) are shown. 




FIG. 10. Total ^Li(n,n)^Li reaction excitation function from 
[23 |. The solid line is the R-matrix fit with the known 3"'' 
state at 2.25 MeV and the new 0^, 1^ and 2^ states at 2.8, 
3.1, and 3.3 MeV. The dashed line shows the effect of shifting 
the 0+ state down by 200 keV. 



biguous. The contemporary (for 1987) shell model pre- 
dictions were used in Ref. [23] as guidance for the fit. We 
performed our own R-matrix analysis of the n-|-^Li exci- 
tation functions and attempted to incorporate the new 
low-lying 0"^ and 2+ states into the n-|-^Li fit. It appears 
that the low energy n-(-^Li excitation function for elas- 
tic scattering can be reproduced with the O"*", l"*" and 2+ 
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FIG. 11. (Color online) Total ''Li(n,n')'^Li(l/2-) reaction ex- 
citation function from [251] . The dashed blue curve is the R- 
matrix fit with the reduced width parameters of Table |T] The 
2^ from the table overestimates the cross section at around 
1.2 MeV in the c.m.s. The solid red curve is the same R- 
matrix fit, with the elastic component of the 2^ reduced to 
better fit the data. 



states if they are placed at excitation energies above 2.8 
MeV without any modifications to their reduced widths. 
The total cross section for the ^Li(n,n) ''Li(g.s.) reaction 
is shown in Fig. (TU] The solid line is the R-matrix fit 
with the states mentioned above, the background states 
from Table m and the known 3+ at 0.22 MeV. Note that 
if the 0+ is shifted down by as little as 200 keV it would 
appear as relatively narrow peak, which is not observed 
experimentally (dotted line in Fig. llOp . From the con- 
siderations above we can conclude that existence of the 
new low lying 0"*" and 2^ states in *Li does not, in prin- 
ciple, contradict available n-|-^Li elastic scattering exper- 
imental data. However, these states have to be shifted 
up in excitation energy by ~800 keV compared to their 
suggested location in ^B. There is also strong evidence 
against degeneracy of the new state(s) with the S"*" state. 
If such degeneracy exists then the experimental cross sec- 
tion at the maximum of the S"*" peak (0.22 MeV) would be 
higher than can be accounted for by the S"*" state alone. 
Our R-matrix fit shows that this is not the case. 

It is more difficult to reconcile the new states in and 
the available '^Li(n,n')'^Li(l/2-) experimental data. If re- 
duced widths parameters from Table [J for these states 
are used then the ^Li(n,n')^Li(l/2-) cross section is over- 
estimated due to too strong contribution from the 2+ 
state (dotted blue curve in Fig. [TT]) . The elastic re- 
duced width amplitude of the 2+ state has to be reduced 
from 0.276 to <0.1 in order to produce a good fit to the 
'^Li(n,n')'^Li(l/2-) data (red solid curve in Fig. AW 
other parameters for the 2+ and also all parameters for 
the 0"*" and 1"*" do not require any modification. The 
''Be(g.s.)+p spectroscopic factor for the 2+ state is al- 
ready small (4%) in ^B but it appears that it needs fur- 
ther reduction to less than 1% in *Li to reproduce the 



^Li(n,n')^Li(l/2 ) data. We do not have a good expla- 
nation for this situation. 

The excitation energy shift of 800 keV between states 
in mirror nuclei (Thomas-Ehrman shift (26l [27j ) is very 
large. While not unique (for example, the shift between 
the 1/2"'" second excited state in ^^O and ^^Na is 730 
keV [i^, H^) it is generally associated with single parti- 
cle structure, where the valence nucleon is in the s-wave 
state. A large Thomas-Ehrman shift results from dif- 
ferent asymptotic behavior of the valence nucleon wave 
function between bound and unbound states in mirror 
nuclei (Nolen-Schiffer effect [13 )• The O^" and 2+ are 
p-shell states, therefore, a large Thomas-Ehrman shift is 
not expected. Realizing that this is an unusual situa- 
tion, we attempted to reproduce the observed p-|-''Be in- 
elastic scattering excitation function without introducing 
the new resonances in ®B but assuming a direct excita- 
tion mechanism of the ^Be first excited state in p+^Be 
scattering. 

These calculations were performed using the coupled- 
channels approach. The potential of Kim et al. [3l| was 
first considered for the bare part of the ^Be interaction. 
The transition potential for the coupling between the 
ground and the first excited state was generated deform- 
ing the bare potential and assuming that these two states 
of the ^Be nucleus are members of a if = 1/2 rotational 
band with a quadrupole deformation length of (52 = 2 fm. 
Besides the transition potential, this procedure gives rise 
also to reorientation terms, which were also taken into 
account in the calculations. The coupled equations were 
solved to all orders using the computer code FRESCO 
[S^ l- In Fig. [T2k ) we show the excitation function for 
a Oc.m. = 146° as a function of the p-|-''Be cm. energy. 
Clearly, the contribution of the direct mechanism is very 
small in this energy range, suggesting that the magnitude 
of the measured inelastic cross section at these energies 
cannot be explained by a pure direct reaction mechanism. 
We performed a second coupled-channels calculation us- 
ing a potential that produces a resonance at these en- 
ergies. This potential was parametrized using a Woods- 
Saxon shape, with radius R — 2.23 fm (deduced from the 
matter radius of the ''Li nucleus), diffuseness a = 0.65 fm 
and the depth adjusted to produce a resonance around 
Ec.m. — 2 MeV. The calculated inelastic excitation func- 
tion obtained with this potential is given by the solid line 
in Fig. WBd ) . The presence of the resonance produces a 
pronounced maximum about 2 MeV and a significant in- 
crease of the magnitude of the cross section. So, based on 
the coupled-channels analysis we conclude that the high 
inelastic scattering cross section cannot be reproduced 
unless resonance(s) is(are) introduced in the correspond- 
ing energy range. 

Finally, we have to make an important distinction be- 
tween the 0"*" and the 2"'" states. While existence of the 
2"*" state is hard to dismiss, the case for the 0"*" state 
is somewhat weaker. In spite of the fact that without 
this resonance the inelastic cross section at 2.0 MeV due 
to direct excitation of the first excited state in ''Be is 3 




FIG. 12. Inelastic scattering differential cross section from 
the ^Be(p,p') reaction calculated within the coupled-channels 
approach, assuming a direct mechanism. The left panel uses 
a p+'^Be potential which does not contain resonances within 
this energy interval. The right panel shows the result of the 
calculation using a potential that contains a resonance. 



times smaller than the experimental value, one should 
be careful making the final call based on such evidence. 
Further investigation is warranted. Specifically, accurate 
measurement of the p-|-''Be excitation function of inelas- 
tic scattering in the energy range from 0.7 to 2.0 MeV 
and in a broad angular range should provide a definitive 
answer on the existence of the 0+. At this point we can 
only regard this state as tentative. 



V. THE CONTINUUM SHELL MODEL 
ANALYSIS OF THE P+^BE DATA. 

The time dependent continuum shell model Q was 
used as an alternative and more microscopically con- 
strained way to analyze the p-|-^Be data. This model 
extends the traditional shell model into the domain of 
reaction physics. It incorporates the many-body dynam- 
ics with all essential structure and reaction components, 
and allows one to predict the reaction observables. Some 
features, such as the angular dependence of cross sec- 
tions and interference between resonances are particu- 
larly sensitive to the many-body structure. The TDCSM 
is built upon one of the well-established Hamiltonians of 
the traditional shell model coupled to reaction channels, 
where a Woods-Saxon shaped potential is taken from a 
global Woods-Saxon parametrization [33]. This theoreti- 
cal treatment of using the WBP shell model Hamilto- 
nian |33| is reported in Ref. The WBP Hamiltonian 
was selected because unlike most interactions it results in 
low-lying 1^, 0^, and 2^ states in ^B, at excitation ener- 
gies below 3 MeV. To consider a full spectrum of possible 
Hamiltonians in this investigation, in addition to WBP, 
we use PWT ^] and CKI [14,] shell model interactions. 
The comparison of the experimental spectroscopic factors 
of the positive parity states in ®B to the predictions of 
the shell model with different interactions is given in Ta- 
ble |lTl The experimental spectroscopic factors were calcu- 
lated as the ratio between the partial width and the single 
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FIG. 13. (Color online) Inelastic scattering differential cross 
sections for ^Be(p,p') reaction obtained with the TDCSM that 
uses three different Hamiltonians is compared to the experi- 
mental data. The panels (a), (b), and (c) correspond to WBP, 
PWT, and CKI interactions. 



particle width calculated using a Woods-Saxon potential 
with a global Woods- Saxon parametrization [3^. In Ta- 
ble|lT]it can be seen that all three residual interactions are 
in good agreement with the experimental spectroscopic 
factors for the Ij'" and 3^ resonances. All three interac- 
tions reproduce the inelastic spectroscopic factor for the 
0^ , and both the elastic and inelastic components of the 
ij. The WBP interaction is the only interaction that 
predicts a 2^ that is dominated by an inelastic compo- 
nent, as is seen experimentally, while the PWT and CKI 
interactions both predict a 2"*" state with a similar inelas- 
tically dominated component as the 2'^ . 

The best validation of the theoretical model predic- 
tions can be performed if the measured cross section is 
calculated directly from the model. Unfortunately, the 
reaction physics is very sensitive to kinematics and to 
the exact position of levels in the spectrum because of 
the phase space and barrier penetrability. While the tra- 
ditional shell model may, in general, be good in describing 
positions and ordering of states, often its precision is not 
close to what is required by the reaction physics. Thus, it 
is common practice to set the exact reaction kinematics 
based on observation. In our approach all known states 
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TABLE II. Experimental spectroscopic factors of positive par- 
ity states compared to the shell model predictions. 
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and thresholds are adjusted from experimental data and 
we treat the energies of unknown ij, 0^, and 2^ states 
as parameters. In our study we vary these three param- 
eters to best fit the observed cross section. The TDCSM 
provides an effective mechanism to modify the position 
of any state in the Hamiltonian while keeping all struc- 
tural aspects unchanged. This is done by adding to a 
shell model Hamiltonian a factorisable term 6E\a){a\, 
where \a) is the eigenstate to be shifted and SE is the 
energy shift. The corresponding change in the many- 
body propagator is performed exactly with the help of 
Dyson's equation, for details see Ref. jSl- 

In Fig. [T2] the inelastic scattering cross section for 
^Be(p,p') obtained with TDCSM is compared to exper- 
iment. Panels (a), (b), and (c) correspond to calcula- 
tions with WBP, PWT, and CKI interactions, respec- 
tively. The spin and parities of resonances in the energy 
region plotted are marked. The inelastic cross section is 
not sensitive to the S"*" state which is seen in the elas- 
tic scattering cross section. All models predict a similar 
structure of the 3+ state and therefore produce a compa- 
rable elastic cross section which agrees well with experi- 
ment. The elastic ^Be(p,p) cross section with WBP inter- 
action is demonstrated in Ref. Q. Positions of ij, 0^, 
and 2^ resonances, indicated in Fig. 1131 are not known 
a priori; here they are adjusted by visual examination to 
best reproduce the experimental data. The main peak in 
the ^Be(p,p') cross section is due to the 2^ resonance at 
around 2.5 MeV of excitation. It was found that agree- 
ment with the experimental data is good if the Of is 
placed around 2 MeV and the 1^ is moved to higher ex- 
citation energy. (Sensitivity of the inelastic cross section 
to the position of the 1^ state is weak. However, the 1^ 
state at excitation energies below 2.3 MeV would produce 
a peak in the elastic cross section, which is not observed 
experimentally. See discussion in Section iHll ) In the case 



of the WBP interaction. Fig. rTST a'). no position adjust- 
ment was made to the 1^ and 0^'" states and the 2^ is 
only moved down by about 140 keV. The CKI Hamilto- 
nian gives two 2+ excited states at 4.2 and 5.1 MeV of 
excitation. Both of these states have been tried as candi- 
dates for the 2.5 MeV resonance and it was determined 
that the second 5.1 MeV state in the CKI Hamiltonian 
has the correct structure. Our main conclusion from the 
calculations shown in Fig. [13] is that the CKI interaction 
appears to be best in reproducing the cross section. The 
states obtained with the CKI appear to have structure 
which agrees with the observed interference and angular 
dependence features. In particular, only the CKI inter- 
action is able to reproduce the observed increase in the 
cross section at 2.5 MeV for higher angle (Fig. [T^ c)). 
The height of the resonance peak at 2.5 MeV is the pri- 
mary difference between theory and observation. We at- 
tribute this difference partially to the 1" state, which 
was not included in the shell model analysis (only p-shell 
states were considered) and also to the somewhat differ- 
ent ratio between the elastic and inelastic partial widths 
for the 2 J. 

The amplitudes from the final R-matrix fit with the 
CKI interaction are summarized in Tab. IIIII (Resonance 
reduced width parameters from the R-matrix fit were re- 
coupled from the LS to the jj coupling scheme for di- 
rect comparison with TDCSM amplitudes). It should be 
noted that the choice of channel radius in the R-matrix 
calculations will have the affect of scaling the reduced 
width parameters, thus one should not directly compare 
absolute values, but rather sign and relative values of the 
reduced widths. We note that the fit only slightly modi- 
fies the amplitudes for the 3^ and 2^ states leaving the 
general features of CKI unchanged. However, there is 
a significant difference between experimental excitation 
energies and CKI predictions for the 2^ state. The 0^ 
also has a significant shift between the experimental and 
CKI predicted energy. 



VI. COMPARISON OF THE AB INITIO MODELS 
TO THE NEW EXPERIMENTAL DATA. 

The first ab initio calculations for A=8 isotopes were 
performed in 1998 by P. Navratil and B. Barrett [l^ us- 
ing the large-basis no-core shell model approach. Except 
for the broad l"*" state at 3.2 MeV (which was under- 
bound by 2 MeV) the known states have been well re- 
produced by the calculations. In addition to the known 
states the new 0+ and 2+ states have been predicted 
at 4.5 and 5.0 MeV respectively. While these states had 
rather high excitation energy they were still the third and 
the fourth excited states in the calculations, below the 
second 1+ state (predicted at 5.2 MeV). Refined NCSM 
calculations with CD-Bonn 2000 potential [s^ produced 
the O"*" as the second excited state at 2.23 MeV in ®B and 
2.48 MeV in ^Li. The 2+ state was produced at 3.8 MeV. 
This prediction is remarkably close to the experimental 



11 



TABLE III. Decay amplitudes from the final R-matrix fit and from the CKI based TDCSM. The first column denotes the spin 
and parity of the resonance, following are excitation energy and four amplitudes from the R-matrix fit. Excitation energy and 
amplitudes for states in *B from the CKI Hamiltonian are listed in the remaining five columns. 
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0.53 
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result of this work with respect to the 0+ state, suggested 
at 1.9 MeV in ^B. Another interesting prediction is that 
the excitation energy of this state should be shifted up 
in ^Li by 250 keV. Much smaller shifts were predicted for 
the known 1+ and 3+ states (50 and 40 keV respectively). 

Similarly to the NCSM, the new low-lying positive par- 
ity states were predicted by the Green Function Monte 
Carlo method [la]- The 0+ state was predicted as the 
second excited state with excitation energy of 1.91(29) 
MeV (the results were given only for ^Li). The GFMC 
results for the 2+ state are not available, nevertheless, the 
variational Monte Carlo calculations produce this state 
at 3.86(18) MeV However, the more recent version 
of CFMC calculations ^ with the AV18/IL2 Hamiho- 
nian puts the 0+ and 2+ at excitation energies of 3.6 and 
5.3 MeV, respectively. (This version of the GFMC calcu- 
lation also puts the known 1+ state at 4.7 MeV instead 
of 3.2 MeV.) 

Clearly, there is no unified picture for the level struc- 
ture of the ®B - ^Li isotopes from the available array of 
ab initio calculations. However, all of them produce a 0+ 
as either the second or third excited state, always below 
the known 1^ state (experimentally found at 3.2 MeV in 
*Li). Our experimental result confirms this prediction. 
The 2+ state is generally found at higher excitation en- 
ergy in ab initio calculations than observed in this work 
for ^B. 

With the development of the ab initio NCSM/RGM 
approach Q, the ab initio phase shifts can now be com- 
pared directly to the experimental phase shifts extracted 
from the R-matrix analysis of the experimental data. The 
''Be-l-p diagonal phase shifts as well as ''Be(p, p') excita- 
tion function have been calculated in Ref. [l3|- Fig. 
1141 shows the ab initio negative parity 1~ and 2~ phase 
shifts from [3] as red and black dashed curves respec- 
tively. The experimental 1" and 2~ phase shifts from 
the R-matrix best fit are the red and black solid curves 
respectively. The experimental and the theoretical phase 
shifts appear to be different. However, the general trend 
is reproduced. Moreover, the fit is not very sensitive to 
the 1~ phase shift. It is possible to make a good fit with 
X^=0.92 using the 1" phase shift from While it 
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FIG. 14. (Color online) Phase shifts from the R-matrix best 
fit and from Navratil [13] for the negative parity states. R- 
matrix calculations are the solid black and red curves for the 
2~ and 1~ respectively. The calculated ab initio phase shifts 
from 17] are the dashed black and red curves for the 2~ and 
1" respectively. Note that the sensitivity of the fit to the 1~ 
phase shift is very weak as discussed in Chapter I VII 



is marginally worse than the best fit (0.89) it does 
not differ significantly visually and the parameters for all 
other states remain within uncertainty quoted in Table ID 
Therefore, the sensitivity of the fit to the 1~ phase shift 
is small and seemingly large difference shown in Fig. [14] 
should not be considered as significant. The fit is more 
sensitive to the 2" phase shift. It does not reach the max- 
imum value of 12 degrees predicted by the NCSM/RGM. 
If the 2" phase shift from ^] is used in the R-matrix fit 
then the of the best fit is 2.2 and the fit is visually 
worse. This is shown in Fig. [15] by the dashed green 
curve. 

The shape of the 3+ diagonal phase shift is determined 
by the 3+ state at 2.29 MeV, which only has contribu- 
tion from channel spin S=2 and no inelastic component. 
NCSM/RGM overestimates the excitation energy of this 
state by ~1 MeV as can be seen in Fig. [16] otherwise 
the phase shift would be similar to the experimental one. 
This was verified by shifting the experimental excitation 
energy of the 3+ state to the NCSM/RGM result in the 



12 





140 




120 




100 






(mb/ 


80 


a 


60 


"B 




-a 


40 








40 




35 




30 


in 










20 


a 




T3 


15 


~B 


T3 


10 




5 








1 1 


1 1 


1 1 1 
. ^Be(p,p) 13^ 


^ \ 

-^Be(p,p') 146° Jll - 

1 1 


\ \ 

-^B^^^138^^^^ - 

1 1 


\ \ \ 

1 1 1 



1.5 



2.5 1.5 



/leV) 



2.5 1.5 2 2.5 3 

Eexc(MeV) 



3.5 



FIG. 15. (Color online) The excitation functions for the ^Be+p elastic and inelastic scattering. The best fit is the solid red 
curve. The dashed green line is the R-matrix fit with a 2~ phase shift matching the work of T?|. The blue dotted curve is the 
fit with the l""" phase shifts matching the phase shifts from |17l] . 
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FIG. 16. (Color online) Phase shifts from our R-matrix best 
fit for the (solid black curve) and calculations of Navratil 
(dashed dotted red line). 



R-matrix calculations. 

Comparison of the 1+ phase shifts from the best fit 



and from the [17[ is shown in Fig. [TTl The phase shifts 
appear to be very dissimilar. The 1^ state shows up pre- 
dominantly in the S=l channel in the best fit (and also in 
the shell model calculations and in tI) and it is located 



at lower energy than predicted in [T3|. This makes the 
contribution from the inelastic channel negligible and the 
S=l phase shift goes through 90 degrees, unlike in [l7| . 
More important difference is that the best fit S=2 phase 
shift barely shows any sign of the 1^ state while the ah 
initio S=2 and S=l phase shifts have about equal contri- 
bution from the \\ state. In principle, a good fit can be 
achieved with the phase shifts similar to those calculated 
in Ref. [3, [l3 (provided that the first excited state is 
shifted down by about 300 keV from where it appears in 
the ah initio calculations). All excitation functions are 
reproduced with this solution, except for the high and 
low angles in low-energy data from Q. This is shown 
in Fig. [S] We believe this discrepancy is significant and 
favors the CKI shell model prediction over the ah initio 
calculations for the structure of the \\ state in *B. 

The 1^ state is responsible for the behavior of the 1+ 
phase shifts above 1 MeV. Comparing the best fit and ah 
initio phase shifts one can notice that at higher energies 
the best fit S=l phase shift is similar to the S=2 ah initio 
phase shift and the best fit S=2 phase shift is similar to 
the S=l ah initio phase shift. We have produced another 
fit using the ah initio 1+ phase shifts from JJTj and varied 
parameters for all other states. This fit has the at 1.2. 
The fit to the inelastic data is visually identical. Qual- 
ity of the fit to the elastic scattering data is somewhat 
worse (see Fig. [T5)) . Observable parameters for all other 
positive parity states were still within the uncertainties 
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FIG. 17. (Color online) The 1^ phase shifts from the best 
fit and from [13]. Best fit R-Matrix phase shifts for channel 
spins 1 and 2 are black and red solid curves respectively. The 
black and red crosses are S—l and S=2 1^ phase shifts from 




-50 











'Be(p,p) 2* 1=1 






'Be(p,p) 2' 1=2 






■■■'Navratil2*l=1 




- 


■■■'Navratil2* 1=2 






















< 1 < < . < 1 . < < . 1 < . < < 1 < . . < 1 < < . 





0.5 1 1.5 2 2.5 3 3.5 4 4.5 

Ecm (MeV) 



FIG. 19. (Color online) Comparison of the experimen- 
tal best fit 2"*' phase shifts to the calculated phase shifts of 
R-Matrix calculations for channel spins 1 and 2 are red 
dashed-dotted and black dashed-dotted curves respectively, 
while those of [13] are solid red and black curves respectively. 
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FIG. 18. (Color online) The 0^ phase shifts from the best 
fit (solid black curve) and from [13| (red crosses). 



quoted in Table U] Generally, while the best fit 1+ phase 
shifts look very different, we find that the data can be re- 
produced reasonably well with the l"*" phase shifts from 
[17| . This ambiguity can be resolved if wide range of an- 
gles and energy is measured with high accuracy and/or 
experiment is performed with the polarized target. 

The 0+ phase shift is defined by the 0+ resonance at 1.9 
MeV. Fig. [IS] shows the 0"*" phase shifts from the best 
fit and the ab initio calculations [l7j . The two phase 
shifts are very similar, indicating that the structure of 
the 0+ state is well reproduced in [l^l- In order to have 
a perfect match between our phase shift and that from 
17| it is necessary to increase the total width of the state 
to ~1 MeV. A fit to the experimental data with the O"*" 
phase shift from [T^] produces a value of 0.97 and 
is almost indistinguishable visually. This is because the 



stronger 0^ is compensated by the modifications to the 
negative parity shifts and the parameters for the positive 
parity states remain almost unchanged. This is why the 
0"*" state has large uncertainly for its widths (see Table 

The only obvious difference between the experimental 
data and the resuhs of NCSM/RGM calculations 
is related to the 2+ phase shift. The 2^ state is pre- 
dicted at 5 MeV by the NCSM/RGM calculations and 
has about equal contribution from both channel spins. 
The R-matrix fit to the observed 2+ state favors channel 
spin S=l. The dominant decay mode for the experimen- 
tal 2^ state is into the first excited state of ^Be(l/2^). 
This produces the characteristic shape of the S=l 2+ 
phase shift that is very different from the ab initio phase 
shifts as seen in Fig. UM It appears that the 2+ state 
predicted in Ref. [3, [l7| and the observed 2+ are two 
different states. We can speculate that the situation here 
may be similar to the predictions of the conventional shell 
model CKI Hamiltonian, that produces two 2+ states at 
4.2 and 5.1 MeV and only the latter has the correct struc- 
ture (see also discussion in Section |Vl. It is possible that 
the lowest 2+ state predicted by the NCSM/RGM calcu- 
lations is not the one observed in this experiment. 



VII. CONCLUSION 

The excitation function for p-|-''Be elastic and inelas- 
tic scattering was measured in the energy range of 1.6 
- 3.4 MeV and angular range of 132 - 148 degrees. An 
R-matrix analysis of the excitation functions indicates 
that new low-lying states in *B have to be introduced in 
order to explain the large inelastic cross section with a 
well-defined peak at 2.5 MeV. These new states are sug- 
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gested to be the 0+ at 1.9 MeV, 2+ at 2.5 MeV, and 1+ 
at 3.3 MeV with width 530 keV, 270 keV, and 3.2 MeV 
respectively. Evidence for the 2+ state at 2.5 MeV is 
reliable. The 1+ is needed to reproduce the high-energy 
inelastic cross section and is seen in the mirror nucleus 
^Li. The 0+ at 1.9 MeV can be considered as tentative 
and more accurate measurements are needed, especially 
at the lower excitation energy region. However, uncer- 
tainty is not related to the spin- parity assignment. If 
there is a state at 1.9 MeV then it has to be the 0+, as 
any other spin-parity assignment does not allow fits to 
the elastic and inelastic scattering simultaneously. The 
uncertainty is related to the possibility of explaining the 
observed enhancement of the inelastic scattering cross 
section at energies below 2.3 MeV by direct excitation 
of the first excited state in ^Be-|-p scattering. Coupled- 
channels calculations of the ^Be(p,p') inelastic scatter- 
ing cross section assuming a direct mechanism have been 
performed and it was shown that without resonance(s) 
the cross section is significantly smaller than observed 
experimentally. Nevertheless, taking into account uncer- 
tainties of the coupled-channel calculations we consider 
the 0+ state as tentative. 

Analysis of the available experimental data on the mir- 
ror nucleus, *Li, indicates that it is unlikely that the new 
states can be found at excitation energies below 2.8 MeV 



in ^Li. Therefore, the excitation energy shift for these 
new states between the two mirror nuclei is ^800 keV. 
This is a factor of two or three larger than the typical 
Thomas-Ehrman shift in p-shell nuclei. We cannot offer 
an explanation for this phenomenon. 

The Time Dependent Continuum Shell Model provided 
important guidance in the R-matrix analysis of the exper- 
imental data. Due to the presence of many broad over- 
lapping resonances the number of free parameters is large 
and a blind R-matrix fit is ambiguous. Therefore, theo- 
retical constraints become very important for extracting 
more reliable results from the R-matrix fit. This is a 
typical situation for the resonance scattering with exotic 
nuclei and development of the theoretical tools for this 
problem is an important step forward. 

The recent development of ab initio NCSM/RGM cal- 
culations open up an exciting opportunity for a sensitive 
tests of ab initio models. We performed a detailed com- 
parison of the diagonal p-|-^Be phase shifts calculated in 
Ref. d, [l3l to the experimental data. Overall, with the 
exception of the 2+ phase shift the results are encourag- 
ing and the predictions of the ab initio model are close 
to the R-matrix best fit. 

The work was supported by National Science Founda- 
tion under grant number PHY-456463 and by the U.S. 
Department of Energy grant DE-FG02-92ER40750. 
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